Mon. Not. R. Astron. Soc. 000.m[T7l(2004) Printed 2 February 2008 (MN KTpX style file v2.2) 



Baryon Oscillations and Dark-Energy Constraints from Imaging 
Surveys 



Derek Dolney 1 *, Bhuvnesh Jain 1 *, Masahiro Takada 2 * 

1 Department of Physics and Astronomy, University of Pennsylvania, 209 S 33rd St, Philadelphia, PA 19104-6396, USA 
2 Astronomical Institute, Tohoku University, Sendai, 980-8578, Japan 



2 February 2008 



ABSTRACT 

Baryonic oscillations in the galaxy power spectrum have been studied as a way of probing 
dark-energy models. While most studies have focused on spectroscopic surveys at high red- 
shift, large multi-color imaging surveys have already been planned for the near future. In view 
of this, we study the prospects for measuring baryonic oscillations from angular statistics of 
galaxies binned using photometric redshifts. We use the galaxy bispectrum in addition to the 
power spectrum; this allows us to measure and marginalize over possibly complex galaxy bias 
mechanisms to get robust cosmological constraints. In our parameter estimation we allow for 
a weakly nonlinear biasing scheme that may evolve with redshift by two bias parameters in 
each of ten redshift bins. We find that a multi-color imaging survey that probes redshifts be- 
yond one can give interesting constraints on dark-energy parameters. In addition, the shape 
of the primordial power spectrum can be measured to better accuracy than with the CMB 
alone. We explore the impact of survey depth, area, and calibration errors in the photometric 
redshifts on dark-energy constraints. 

Key words: cosmological parameters - distance scale - equation of state - large-scale struc- 
ture of Universe 



1 INTRODUCTION 

Recent cosmological observations have provided strong evidence that a dark-energy componen t, such as the cosmological constant, com- 
prises as much as 7 per cent of the total energy density of the universe IPerlmutter et alll999tlRiess et alll99il200ll:lTonrv et alfcooi 
ISDergeletall2 003). Characterizing the nature of the dark energy and its possible evolution has become a central goal of empirical work in 
cosmology. 

Galaxy surveys measure the clustering statistics of galaxies as a function of scale and redshift. The galaxy power spectrum can be 
compared to the CMB power spectrum to constrain the growth of structure. However the amplitude of the galaxy power spectrum depends on 
the biasing of the particular galaxy sample; one therefore needs to exercise care in using the full power spectrum for cosmological inferences 
(e.g. lTegmark et all2004l : |Pope et all2004ll . The shape of the power spectrum has been regarded as more robust to biasing effects. 

The baryon oscillations in the galaxy power spectrum are imprints from acoustic oscillations in the early universe, prior to recombination. 
The same physics produces the dramatic peaks and troughs seen in the CMB power spectrum. Th eir physical scale is set by the sound 
horizon at recombination, which has been determined accurately from CMB data iSpergel et all2003l) . The baryon oscillations in the matter 
and galaxy power spectra are much weaker features because the dark matter which dominates the mass density did not participate in the 
acoustic oscillations. The oscillations are at the level of a few tens of a percent variation about the smoothed power spectrum. With a survey 
of sufficient size, these features can be measured accurately. Since the physical length scale of the baryon oscillations is known from the 
CMB measurement of the sound horizon, a measurement of their apparent size in redshift or angle space leads to a measurement of purely 
geometric quantities: the Hubble parameter and the angular diameter distance, respectively. We will be concerned with the relation between 
physical size and angular diameter distance: r p h ys icai = <1a(z)A8, where <1a is the angular diameter distance and A0 is the measured 
angular scale. This relation can be used for a statistical quantity as well; for the power spectrum it means that a measurement of an angular 
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wavenumber I and its relation to the physical wavenumber k yields <1a{z). We describe in the next section how d,A constrains models of dark 
energy. 

To measure baryon oscillations, many authors have considered galaxy surveys over dif ferent redshift ranges ( Eisenstein 2003l lBlake & Glazebrool 
l2003tlLindeJ2003l : lHu & Haimarl2003tlSeo & Eisensteml2003l : lMatsubara & Szalavl2 003). For spectroscopic redshift surveys, the tangen- 
tial and radial components are considered separately since the latter is subject to redshift distortions. Current redshift surveys can map large 
enough volumes at redshifts well below 0.5. It is a great advantage to probe higher redshifts since additional baryon oscillations can then be 
measured within the linear regime of clustering (the linear regime extends to smaller physical scales at high redshifts). With future redshift 
surveys, such as the proposed KAOS 1 survey, such a measurement would be possible. 

Multi-color imaging surveys are already in progress, e.g. the SDSS 2 , the CFH Legacy 3 survey, and deeper surveys are proposed for 
the future, e.g. PANSTARRS 4 , LSST 5 , SNAP 6 and others. These surveys offer the possibility of photometric redshifts as crude redshift 
estimators. With the SDSS providing a large sample of relatively nearby galaxies, and the Hubble Deep Fields 7 and the GOODS 8 survey 
providing deep samples of galaxies beyond z = 1, many multi-color samples of galaxies have been studied and used to estimate photometric 
redshifts. With good photometry in 4—5 optical filters, it is expected that a statistical accuracy characterized by an rms of o z ~ 0.02 — 0.04 
in 1 + 2 is achievable for galaxies below z = 1. For special samples such as the Large Red Galaxy (LRG) sample of the SDSS, one can do 
significantly better. Similarly it is expected that with more filters and experience with estimating pho to-z's, the accuracy will im prove and 
extend to higher redshifts. This is an area of ongoing study (e.g. IConnollv et alkoollBudavari et alfcoollMobasher et all2004h . 

The accuracy of photometric redshifts determines the bin width in redshift within which angular power spectra can be measured and 
regarded as being independent of neighboring bins (i.e. the overlap in the galaxy distribution bet ween neighboring bins i s small). This is 
important because wide bins would cause the baryon wiggle features to be smeared out. Following lSeo & Eisenstein ^2003) we will assume 
that a z < 0.04 in 1 + z. Note that this is not a very stringent requirement; at z = 1, it means the rms error in the photometric redshift is 
below 0.08. Given a large number N of galaxies with photo-z's, the mean redshift is measured accurately since the error in it is suppressed 
by vN, which can be very small even per redshift bin for surveys of several 100 or 1000 square degrees. However, when the photo-2's are 
estimated in terms of a probability distribution over redshift per galaxy, often the result is bi-modal or worse. Thus there are sometimes large 
errors in the estimated photo-2, and for certain redshift ranges they lead to systematic biases. While calibration with spectroscopic redshifts 
of some fraction of the galaxies can be used to limit such a bias, we will explore the sensitivity of our dark-energy constraints to possible 
biases in the mean bin redshift. 

Our focus will be on the question: can subtle effects in the biasing of galaxies compromise the dark-energy constraints obtained from 
them? We will use the bispectrum, the Fourier transform of the three-point function, as an independent measure of the biasing of galaxies. The 
idea of using the bispectrum in addition to the power spectrum on large scales to cons train both the amplitude of the mass power spectrum 
and bias has been suggested and implemented iFrv 1999; Frieman & Gaztanasa 1999; Feldman et al. 2001; Verde et al. 2002). In this study 
we will examine whether it can help constrain a possibly redshift dependent bias that could mimic effects of the cosmological model we are 
trying to constrain. We will also compute the gain in signal to noise that the bispectrum provides. We will use perturbation theory, valid in 
the weakly nonlinear regime of clustering. 



2 THEORY 

2.1 Dark-Energy Cosmology 

The expansion history of the universe is given by the scale factor a(t) in a homogeneous and isotropic universe. For cosmological purposes, 
dark energy is associated with a density field, possibly weakly evolving over cosmic time, that affects the expansion rate of the universe, 
H(a). The expansion rate at late cosmic times depends on contributions from the density of non-relativistic matter, fi m (cold dark matter 
plus baryons), and the dark-energy density, Qdc- Our notation is such that fi m and Q<ie denote the present day values in units of the critical 
density, 3H% /(8nG). The dark -energy formalism introduces a new, time-dependent density described by the dark-energy equation of state: 

, \ Pdc 1 dlnpde 1 

w(a) := = -g-Ti L (1) 

Pde o a In a 

This density component is allowed to have an arbitrary time dependence; it must be constrained by theory and observation. Due to lack of 
compelling theoretical models for the dark energy, it is typical to introduce a parameterization for w(a) to be determined empirically. We 
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follow lLinde3l2003l) and parameterize the equation of state as 

w(a) = wo + w a (l — a). (2) 

Note that setting w = — 1, i.e., wo — — 1 and w a = 0, corresponds to the equation of state for a cosmological constant. 
We assume a flat universe throughout this paper. The expansion rate for a fiat universe with dark energy is given by 

if 2 (a) = J ff ( ?|fi m a- 3 + n do e- 3/ i adlna '[ 1+ra(a ' ) ]] . (3) 

The modification in the expansion rate of the universe by a dynamically evolving dark-energy component, as opposed to a cosmological 
constant for example, alters distance measurements on cosmological scales. The (comoving) distance, x( a )> from an observer at redshift zero 
to a source at scale factor a — 1/(1 + z) is given by: 

Hence, given a cosmological standard ruler, such as the length scale associated with one of the peaks of baryon-induced oscillations in the 
power spectrum at a given redshift, one may obtain dark-energy constraints via constraints on the distance function. 

Dark energy that has negative pressure, i.e., pdc + 3pdc < as observed today, leads to repulsive gravity and therefore does not cluster 
significantly. True, we are interested in cluster ing on la rge scales in this pa per, yet spa tial fluctuations in the dark energy are not expected to 
be significant below the present horizon scale ICaldwell et alll999l:lMa et all200rilH32002l) . which we do not approach. We thus treat dark 
energy as a function only of epoch. 

However, since the expansion rate is altered, time-dependent dark energy does modify the redshift evolution of mass clustering. In linear 
theory, all Fourier modes of the mass density perturbation, S := 8p m /pm (relative to the ), grow independently and at the same 

rate, so Sk(a) =: D(a)Sk(a = 1) defines the growth factor, D(a). The growth factor is obtained by solving the linearized perturbation 
equation: 8k + 2H5k = 47rGp m <5fc, where the dot denotes derivative with respect to physical time. This equation can be manipulated to 
obtain a differential equa tion for the growth factor or, alternatively , the growth rate relative to that for a flat, matter-dominated universe, 
g(a) := D(a)/a JWang & SteinhardJl998|.lLinder & Jenkinsl2()03l) : 

?f + 5 P- 3w(a)n de (a)] = § [w(a) - 1] Sl do (a)4r. (5) 
atr Z a da A er 

The dark-energy density parameter at epoch a is denoted Side (a), and is obtained from the relations 

" dc(a) =iTTR' and (6) 

X(a) = e-3.fi dW^(aO (7) 

1 i z m 

_We restrict our attention to the quasi-linear regime and expand the matter density contrast perturbatively as te.g.. lBernardeau et alJ 

120 



S = D8i +D 2 S 2 +0(S 3 ). (8) 
We suppose a biasing prescription to obtain the galaxy density contrast from the matter density: 

S s := ^ (9) 

= M+ h 2 8 2 +0(S' A ). (10) 

We allow for primordial tilt and run in the line ar power spectrum b y introducing the spectral index, n s , and run, a := din n s /d In k, 
parameters defined at ko = 0.05/Mpc, as in Kosowskv & Turner 1 1995): 

/ u \ «s + |aln(fe/fe ) 

P h (k)=A - T 2 (k). (11) 



For T(k), we use the transfer function o f lEisenstein&HuNl998h that includes baryon oscillations. We determine A by CMB normalization 
using the relation 

~ 4 r a. i 2 „ / kt \- J 

(12) 



A 4 A 2 ( H ° 



2tt 2 



That is, the normalization is determined by normalizing the density perturbation amplitude at fco = 0.05/Mpc from the primordial curvature 
fluctuation contrast, 8( . The epoch 2* should be representative of the matter-dominated regime. The normalization is quite insensitive to this 
choice. We take 2, = 1088. We will employ the CMB scalar fluctuation amplitude parameter, A s , whose relation to is 

A s = (1.845 c x 10 4 ) 2 . (13) 
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z 



Figure 1. Galaxy redshift distribution and binning scheme employed in our analysis. The j/-axis gives the number density of galaxies per redshift bin of width 
Az = 0.1. The numbers in the shaded regions indicate the unused galaxy number density in the low-z and high-z regions. The square symbols show the 
signal-to-noise for the power spectrum and the crosses for the bispectrum (see right side of the j/-axis) 



2.2 Galaxy Statistics 

We calculate the projected galaxy power spectrum using the Limber approximation: 



P(l) 



d XPg(x) 



1 



where a — 1/(1 + z) is the scale factor, 

P g (k,a) = blD 2 (a)P L (k), 

and D(a) is growth factor of the previous section. 

For the distribution of sample galaxies, p g (x)> we use the source distribution of lrluterej i2002t) : 



-z/z 



(14) 



(15) 



(16) 



with zq — 0.5. We separate this distribution into 10 bins over the redshift range from z = 0.3 to 1.3. All bins have equal width (0.1 in 
redshift), as indicated in Fig.Q We project the galaxy distribution for each bin of galaxies separately using Eg , 1141 A more rigorous treatment 
would convolve the galaxy distribution with an appropriate selection function for a galaxy catalog with photo-z information. Our approach 
is a bit more crude, but we attempt to account for redshift uncertainties by adding a nuisance parameter for the mean redshift of each bin to 
our Fisher matrix, as described in the next section. 
Since galaxy number is conserved, we may write 



Ps(x) d X = Pg{z)Az. 
Finally, we rewrite Eq. 1141 as 
1 



P\l) 



(n|)3 



Ps(2) dx" 



X 2 



-, a, 
X 



(17) 



(18) 
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Parameter 


Fiducial Value 


Description 


Calculated As 




0.023 


Baryon Physical Density 




w d 


0.112 


Dark-Matter Physical Density 




fide 


0.69 


Dark-Energy Density Parameter (Eq.l3l 




W 


-1 


Equation of State Parameter (Eq.l2l 




w a 





Equation of State Parameter (Eq.l2l 




A s 


0.82 


Scalar Fluctuation Amplitude at k = 0.05/Mpc 






0.979 


Primordial Spectral Index at k = 0.05 /Mpc (Eq.I1 11 




ex 





Primordial Run (= d\nn s /dlnk) at k = 0.05/Mpc 




hi 

u 1 


0.998 


First Order Galaxy Bias Factor (Eq 




hi 





Second Order Galaxy Bias Factor fEq.HOl 




T 


0.143 


Optical Depth 






0.31 


Matter Density 




n h 


0.053 


Baryon Density 


ui h /h 2 


fitot 


1 


Total Density (Assume Flat Cosmology) 




h 


0.66 


Current Hubble Parameter in Units of 100 km / s / Mpc 


\/K +W d )/(1 - fide) 


°"8 


0.88 


Galaxy-Scale Fluctuation Amplitude 





Table 1. The fiducial cosmological model used throughout this paper. 



where we have introduced the index i to label the different redshift bins. 
We calculate the full-sky bispectrum from the flat-sky bispectrum as 



h h h 




(2l 1 + l)(2l 2 + l)(2l 3 + l) 

47T 



B(li,h,h), 



with the flat-sky bispectrum likewise given by a Limber-style projection: 

1 3 , 



(n|)3 



dz 
Ps(z) d£. 



■ I h h h 
—tif, | a 

XXX 



(19) 



(20) 



To 0(8 i ), the galaxy bispectrum is given by the sum of two terms: 

B s (ki, k2, k^,a) = biD 4 (a) [2i*2(fci, k2)P(k\)P(k2) + cyclic ki permutations] 
+ -61&2D 4 (a) [2P(fci)P(fe) + cyclic ki permutations] , 



where 



F 2 (fei,fca) 



5 1 fei • fea 

7 + 2 kik 2 



ki fa 
ki ki 



+ 



2 (kx-ki 



7 V kik2 



(21) 



(22) 



is the second-order perturbation theory kernel. Note that the two terms have different triangle configuration dependence, and different depen- 
dence on bi and 62. 



2.3 Covariance and Fisher Matrix Formalism 



The power spectrum covariance for redshift bin i ignores the contribution from the trispectrum, which should be small in the linear regime: 

, i 2 



/sky 21 + 1 



(23) 



Here / s k y is the fraction of sky coverage for the survey, for which we use the value 0.1. n 1 ^ is the galaxy density for redshift bin i. In our 
galaxy distribution (Fig.0 it is 2.4 galaxies per arcmin 2 for the lowest redshift bin, and 5.4 galaxies per arcmin 2 at the peak of the galaxy 
distribution. 

The bispectrum covariance ignores contributions from 3-, 4-, and 6-point functions, which are also expected to be small: 



C, 



m L 



P(h) + — 



p (h) + ±- 



(24) 

We assume that the likelihood function for the galaxy power spectrum and galaxy bispectrum is Gaussian, and employ a Fisher matrix 
analysis to approximate the likelihood near the fiducial cosmological model defined by the parameters in TableQ 
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Figure 2. Angular separation and wavenumber of the first five peaks in the galaxy power spectrum induced by baryon oscillations as a function of galaxy 
redshift. The dotted lines correspond to our cutoff values 2 max and Z m j n . The 1(8) relation is approximate — we use I = 2n/8. 



The Fisher matrix elements are given by 

*5. = £ E f^r 1 ^ ( 25 ) 

for the power spectrum, and 

F «P-1^ -^—^hhh^hh) —q— (26) 

' imin^il^i2^i3SSimax 

for the bispectrum. The derivatives are evaluated at the fiducial model. 

The large-scale limit, Z m i n , is introduced because the Limber approximation fails when probing scales larger than the line-of-sight width 
of the given bin of galaxies. So we let l m i n be the wavenumber that corresponds to the width of the bin. Since the information from baryon 
oscillations is contained at significantly larger I in any case, this choice has little impact on our results. 

The small-scale limit, Z max , is meant to represent the scale beyond w hich perturbation theory is unacceptable. Nonlinear effects tend 
to smooth out the high-Z baryon oscillations in any case (Meiksin, White & Peacockl ll999h . We define this (redshift-dependent) scale as 
imax = fc max /x(z), where fc max is defined by 

' aX d 3 fcP(fc) = 1. (27) 

The values of Z max are plotted as the upper dotted line in Fig. [2] We will test below the sensitivity of our results to the choice of Z max by 
repeating our calculations for Z max /2. 

We neglect the cross covariance between the power spectrum and bispectrum, which is proportional to the 5-point function, so that the 
Fisher information matrix for the power spectrum and bispectrum combined is simply the sum 

Fc/3 ~ Fa/3 + Fa/3 ■ (28) 
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Finally, for a pair of parameters p a and pp, we calculate \ 2 as 



2 



,P/3 



Ap a 

AjJ/3 



(F- 
(F- 



' /3a 



(F- 
(F- 



'a/3 
'/3/3 



Ap,3 



(29) 



and identify the one-sigma contour as that contour defined by \ — 2.3. 

For a flux-limited sample of galaxies, one should allow for a dependence of the galaxy bias on redshift. Hence we include an independent 
b\ and b\ for each redshift bin in our Fisher matrix. 

We allow for biases in photometric redshifts by including the mean redshift of each galaxy bin as a nuisance parameter in the Fisher 
matrix. We set the rms in the mean redshift to be 0.01 in 1 + z for each redshift bin. This is conservative, since samples of spectroscopic 
redshifts will likely be available for calibration and the residual bias should be significantly smaller. 



3 RESULTS 

3.1 Small Parameter Set 

Since we are most interested in constraining the dark energy and galaxy bias, we first consider a restricted set of cosmological parameters: 
the dark-energy parameters, the overall amplitude normalization, and the bias parameters: {wq, w a , A s , 6i, 62}- We will impose no external 
priors on these parameters, but all other parameters are fixed to the fiducial values in TableQwith the exception of fide, which we marginalize 
over after applying a prior of 0.05. It is important to constrain fi<jc to be able to extract information about the other parameters. Fortunately, 
CMB information can be used to constrain Q^c well. In the next subsection, we use a more complete set of parameters and combine our results 
with information from the CMB, and further allow for a possibly redshift-dependent galaxy bias by giving each redshift bin an independent 
b\ and b\. 

Since the dark-energy parameters affect the angular diameter distances, they shift the position of the peaks in the projected galaxy power 
spectrum and bispectrum. In Fig. [3] we plot the angular galaxy power spectrum for our highest redshift bin: 1.2 sC z 1.3. Notice that, 
in addition to peak shifting, the overall amplitude of the fluctuations is also changed due to the growth factor's dependence on dark energy. 
Because of this, we include those parameters that affect the amplitude of the galaxy statistics in our 5-parameter model. We neglect redshift 
uncertainties, and any redshift dependence of the bias parameters. That is, for this case, a single b\ and 62 parameter are used for all 10 
redshift bins. Also, we plot the galaxy bispectrum signal in Fig.|4| which shows how these parameters modify the configuration dependence 
of the bispectrum. 

The results of our Fisher matrix analysis for the 5-parameter model are presented in Fig. [5] These are two dimensional (marginalized) 
projections of 1-er contours in the 5-parameter space. Based on these contours, the addition of the bispectrum information does not appear to 
strengthen the dark-energy constraints substantially. However, constraints on the other three parameters are all significantly improved by the 
bispectrum information. This is not too surprising, for the galaxy power spectrum is proportional to A s b\ fEq. 1181 . and so the galaxy power 
spectrum alone cannot distinguish A B and 61. The degeneracy is broken by the bispectrum, as its dependence is Aib\ (for 62 = 0, Eq, l20t . 
The degeneracy breaking is most dramatic in the A B — 61 plot of Fig. [5] where the power spectrum or bispectrum information alone results in 
unbounded contours, but each contour comes with a different slopes, so that a tight contour is obtained when the two statistics are combined. 

To lowest order in perturbation theory, the power spectrum does not depend on 62- Dependence on 62 enters at next order, proportional to 
the trispectrum. We have neglected the trispectrum contribution relative to the power spectrum, and so there are no power spectrum contours 
for 62 in Fig. 1181 However our formalism allows for this contribution to be included; with the bispectrum we can actually measure 62 and then 
use the full power spectrum up to second order, which thus allows for a nonlinear, scale-dependent bias. The bispectrum enables constraints 
on 62 because, as can be seen from Eg. 1201 it contains two terms with different triangle configuration dependence. Only the second of the two 
terms depends on 62, so observation of the galaxy bispectrum is useful to constrain this parameter. 



3.2 Large Parameter Set 

We have done the Fisher matrix analysis for a comprehensive set of cosmological parameters: {a>b, u>d, ^de, wo, w a , A s , n s , a, b\,b2, t}. 
Refer to TableQfor a description of each parameter and its fiducial value. We use an independent pair of bias parameters for each redshift 
bin. We also include the mean redshift values of our redshift bins as additional nuisance parameters. We set the rms uncertainty to be 0.01 in 
1 + z for the mean of each redshift bin. 

To get a feel for the sensitivity of the power spectrum and bispectrum signals to each of these parameters, we have plotted in Fig. |6| 
andUJfhe derivatives with respect to each parameter, evaluated at the fiducial model, as a function of /. The derivatives are weighted by a 
covariance factor, 1 / \/C, as this is the weighting that enters the Fisher matrices (cf. Eqs. l25l and l26> . Additionally, for the bispectrum signal, 
a sum over all triangle configurations should grant a proportionately higher sensitivity for larger I values. We attempt to roughly account for 
this by weighting the bispectrum derivatives by a factor of I 2 . Thus, squaring each weighted power spectrum and bispectrum derivative and 
summing is approximately (unmarginalized) l/cr 2 for that parameter. 

With so many free parameters, it is necessary to apply some priors or to combine galaxy data with some complementary cosmological 
information, such as the CMB anisotropy. Galaxy information by itself suffers from strong degeneracies, like u>o-w a in particular. The 
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Figure 3. Effects of dark energy and other parameter variations on the projected galaxy power spectrum (Eg. 1181 . Except for the one labelled exception, we 
plot the projected power spectrum for our redshift bin with z mcan = 1.25. In the lower panel, we divide by a smooth spectrum corresponding to zero baryon 
density. This shows the differences between the models more clearly. Note that the shot noise curve for the lower redshift bins will be more than a factor of 
two larger than what is plotted here for z = 1.25. Even so, it is clear that shot noise is not a significant source of error for the survey parameters we have used. 
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Figure 4. Configuration dependence of the projected bispectrum for various parameters choices as in Fig.0 The upper panel uses equilateral triangles, while 
the lower panel uses isosceles triangles. Note that the effect of the bias parameter 62. which leads to nonlinearity in the biasing, on the bispectrum is at the 
same level as that of 61 . Hence we expect that the bispectrum can be used to constrain this parameter. 
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Figure 5. 1-<t constraints from the Fisher matrix analysis of the 5-parameter cosmological model. The fiducial model is the cosmological constant model in 
TablefTI and is marked with an X in the figure. 



marginalized one-sigma constraints from galaxy data alone are weak: wo constrained to ±1 and w a constrained to ±4. However, combination 
with the complementary CMB information gives interesting constraint ellipses, as shown in Fig.[S| or Table [2] We have not explored other 
possible cosmological probes with which to combine angular galaxy clustering. 

We plot the one -sigma constraints obt ained from our Fisher matrix analysis in Fig. [8] We calculate CMB Fisher matrices in the same 
way as was done for lTakad a & Jairj i2004|). CMB t emperature and polarization power spectra and the cross spectrum are calculated using 
CMBFAST version 4.5 iSeliak & Zaldarri agalll996l) . We assume the experimental specifications of the Planck 143 and 217 GHz channels 
with 65 per cent sky coverage. We also assume a spatially flat universe with no massive neutrinos and no gravity waves. We summarize our 
constraint prospects in Table|2| 

As in Fig.|5| the bispectrum helps primarily in constraining the amplitude and bias parameters. This is important because poor constraints 
on the bi and 62 parameters mean that the dark-energy constraints may not be robust against possibly scale-dependent or nonlinear bias. The 
main difference in the level of accuracy achieved compared to the 5-parameter study is due to the bias parameters, since here we allow for 
two independent parameters in each redshift bin. Constraints on the biasing from other observations, or having some basis to reduce the 
parameter set from the twenty we use here, would improve dark-energy constraints by up to a factor of two. 

Fig. |H illustrates the power of a wide area multi-color survey, in that the full set of parameters gets constrained by the galaxy power 
spectrum and bispectrum, in conjunction with the CMB. For example, the constraints on the shape of the primordial power spectrum and on 
dark-energy parameters are interesting enough that increasing survey area or depth becomes appealing. We discuss these below. 

We performed our analysis for a supergravity (SUGRA) inspired fiducial model, with wq — —0.82 and w a = 0.58. A non-cosmological 
constant fiducial model results in stronger dark-energ y constrai nts as a result of the redshift coverage of galaxy observations at relatively low 
redshift, and CMB observations at high redshift (Lindet l2003h . The comparison in Fig.|5|shows an improvement of over a factor of two in 
constraints on dark-energy parameters, though given Planck level CMB information galaxy clustering does not add much. Table 2 shows 
these results, as well as the use of Type la supernovae constraints at the level of SNAP. 
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Figure 6. Plot indicating the sensitivity of the power spectrum to cosmological parameters as a function of I. The curves have been weighted by the covariance 
1/VC, as in Eq.lBl 



4 DISCUSSION 

We have explored the prospects for obtaining interesting dark-energy constraints from the angular galaxy power spectrum and bispectrum. 
The use of the bispectrum adds some new information, though the signal-to-noise is significantly smaller than that of the power spectrum. 
However, its value lies in making dark-energy constraints robust against biasing features that cannot be measured with the power spectrum 
alone. We have shown that using the configuration dependence of the bispectrum, even a nonlinear, redshift-dependent bias characterized by 
20 parameters in our study, can be measured and marginalized over. 

The dark-energy constraints we obtain are weaker than those expected from Type la Supemovae or weak lensing. An ambitious spec- 
troscopic survey could also do better by using three dimensional information (e.g. Seo & Eisenstein 2003; Linder 2003). However it is 
worth emphasising that our requirements of a multi-color imaging survey place no additional burdens over what is needed for weak lens- 
ing. Surveys such as the CFH Legacy survey are already in progress, and several larger surveys are planned for the future. We have made 
conservative choices of redshift binning and the maximum redshift used, to ensure that photometric redshift uncertainties would be well 
below the requirements. Thus we believe that, at the very least, the angular galaxy statistics would provide a useful independent check on 
dark-energy constraints from imaging surveys. As discussed below, it is possible to improve the constraints to levels competitive with the 
methods mentioned above by increasing survey area, depth and quality of photometric redshifts. 

High redshift (z ~ 1) galaxies are important to obtain the parameter constraints at the level advertised here. This is true for several 
reasons. The linear regime extends to smaller scales at higher redshifts because nonlinear evolution progresses to larger scales over time. 
This is evident from Fig|5| which demonstrates that, at higher redshifts, a greater number of the baryon-induced peaks lie in the linear region. 
The three dimensional nonlinear scale, fc max , is a factor of two smaller at z — 1.3 than at z = 0.3. Further, the two-dimensional nonlinear 
scale, / max , depends not only on fc max , but also on the angular diameter distance, which for z = 1.3 is a factor of 3 larger than at z = 0.3. 
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Figure 7. Plot indicating the sensitivity of the bispectrum to cosmological parameters as a function of I. The curves have been weighted by the covariance 
1 / vC, as in Eg. 1261 and by a factor of I 2 as an estimate of the number of triangle configurations. 



For ongoing and future deep imaging surveys the number density of galaxies in these higher redshifts bins is larger. This is helpful in that 
we can select special classes of galaxies with useful properties, such as the luminous red galaxies studied from the Sloan Digital Sky Survey. 
Even if they are a small fraction of the galaxy population, the shot noise contribution is not likely to be a limiting factor as discussed above 
and by Seo & Eisenstein (2003). 

Since we wish to leam about the time dependence of the dark-energy equation of state, it is important to sample our statistics over a wide 
enough range in redshift. In particular, for the equation of state we have used, which is linear in the scale factor, the degeneracy between Wo 
and w a is broken by sampling galaxies over several redshift bins across a range of redshifts. This is illustrated in Fig. llOl and Fig. II II where 
the marginalized 1-cr constraints for the dark-energy parameters are plotted as a function of survey redshift depth. The constraints from a 
survey that extends to z ~ 2 are twice as good as one that goes to z ~ 0.5. 

The right panel of Fig. llOl shows the dependence of dark-energy constraints on survey area. While the errors do not scale simply with 
/ a ky due to the use of CMB priors, there is a significant improvement with every factor of two in sky coverage. Using Fig.fTol our fiducial 
values can be scaled to a range of survey depths and areas. 

The ability to use the deep imaging capability of future surveys to go well beyond z ~ 1 will depend in part on the accuracy of 
the photometric redshifts. With low enough scatter in these and even a very small sub-sample of spectroscopic redshifts for calibration, 
the methods we have explored can provide much better constraints on dark-energy models, especially models with a larger component of 
the dark-energy density at higher redshifts. As discussed by Seo & Eisenstein (2003), the requirements on photometric redshift scatter are 
not severe: an rms of 0.04 in 1 + z is adequate. Here we have further checked that residual biases in the mean redshift of each bin do 
not significantly degrade dark-energy constraints (see Fig. ll2i . Hence it is not unreasonable to hope that the innermost contour of Fig. llQI 
corresponding to the highest redshift bin at z — 2.25, will be achievable for a high quality survey. In addition, a smaller rms in the photometric 
redshifts will mean that we can use effectively smaller redshift bins to improve the constraints. A rigorous analysis of redshift binning would 
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Figure 8. 1-er constraints from the Fisher matrix analysis of the 11 -parameter cosmological model. We have combined our galaxy power spectrum and 
bispectrum Fisher matrices with a CMB Fisher matrix calculated using projected Planck covariances. The fiducial model is the cosmological constant model 
in TablelTI and is marked with an X in the figure. 
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Table 2. Summary of marginalized 1-c prospects. Except where indicated otherwise, P g and B g include galaxies to redshift 1.3. The bias parameters listed 
here are those for the redshift bin with 0.9 $C z $C 1.0. 
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Figure 9. i-cr constraints from the Fisher matrix analysis of the 11 -parameter cosmological model. The A-model is compared to the SUGRA model in the 
two panels. Note the difference in the range plotted for both parameters — the SUGRA model is constrained much better, though the galaxy statistics add less 
information to the CMB for this model. 
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Figure 10. The two panels show the dependence of dark-energy constraints on survey depth and area. The left panel shows constraints in the w a -WQ plane for 
different values of the maximum redshift bin used for the analysis. The right panel shows constraints for different survey areas, given by the values of the sky 
fraction, / sky . 




Figure 11. Dark-energy parameter constraints as a function of the depth of the galaxy survey. These use the galaxy power spectrum and bispectrum combined 
with CMB data at the level of Planck. The curves are normalized to the constraints expected for a survey with nominal depth of z = 1.25. 



© 2004 RAS, MNRAS 0Q0.fT1[T7l 



1 6 Derek Dolney et al. 



o Dependence 





Figure 12. Dependence of the wo constraint on the uncertainty in the redshift bin centers. Plotted is a Wo normalized to the fiducial value (rms error of 0.01 
in 1 + z). These use the galaxy power spectrum and bispectrum combined with CMB data at the level of Planck. 



require the use of cross-spectra and their covariances; alternatively, the galaxy distribution could be treated as three dimensional, with large 
errors in the redshift coordinate (Seo & Eisenstein 2003). 

Besides the baryon oscillation features, dark-energy parameters affect the overall shape and amplitude of the angular power spectrum 
and bispec trum. To sep arate the effects of the oscillations, we calculated galaxy Fisher matrices using the "no-wiggles" transfer function of 
lEisenstein & Hul ll998l) . If we use just the galaxy clustering information, the constraints for the no-wiggles case are about twice as large. 
With CMB priors, the effect is smaller as the CMB enables the amplitude of the spectra to also provide inform ation. We do still have a 
standard ruler in the no-wiggles case: the peak in the power spectrum corresponding to matter-radiation equality ICoorav et alJl200l[l . The 
physical length scale of this ruler is set at decoupling by uj m , which can be well constrained by the CMB. Hence the broad shape of the 
power spectrum adds information on angular diameter distances. The amplitude depends on the bias parameter and the growth factor, which 
is sensitive to the dark energy. Using both the galaxy power spectrum and bispectrum with the CMB, the scalar amplitude A s , the two bias 
parameters, b\, &2, and the power spectrum amplitude in each redshift bin, are constrained sufficiently well. 

We have tested our choices of fiducial bias parameters and the quasi-linear regime. We repeated our analysis by setting Z max to half 
the value used, and find that it leads to a 6-1 1 per cent degradation in the dark-energy constraints for galaxy plus CMB information. Thus 
our results are not very sensitive to the choice of cutoff of the quasi-linear regime and the neglect of non-Gaussian contributions to the 
covariances. 

Since one of the goals of our formalism is to allow for a nontrivial biasing of the galaxies, we also performed our analysis with b% = 1 
as the fiducial value, thus allowing for a nonlinear term in the biasing. It led to improved constraints on most parameters due to a larger 
impact of adding bispectrum information, e.g. the l-a uncertainty in the dark-energy parameters was lower by 10 per cent. More importantly, 
we find that bi can be measured sufficiently well for a wide range of fiducial values, so that an effectively nonlinear, scale-dependent bias 
can be consistently incorporated in our analysis. This validates our conclusion that bias and dark-energy constraints can be simultaneously 
obtained by using the bispectrum. An alternate approach to galaxy biasing and other properties is to use the halo model, which may also 
enable us to probe smaller scales, as discussed in Hu & Iain 1 2003). The analytical models for the bispectrum need to checked and calibrated 
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with iV-body simulations. A large number of independent realizations would be needed since the covariances in the bispectra need to be 
accurately measured as well. 

ACKNOWLEDGEMENTS 

We thank Gary Bernstein, Daniel Eisenstein, Wayne Hu, Martin White, and Alex Szalay for helpful discussions. We are grateful to Eric 
Linder for discussions and for providing his supernovae Fisher matricies. This work is supported in part by NASA grant NAG5-10924 and 
NSF grant AST03-07297. 

REFERENCES 

Bemardeau, F, Colombi, S., Gaztanaga, E., Scoccimarro, R., 2002, Phys. Rep., 367, 1. 
Blake, C, Glazebrook, K., 2003, ApJ, 594, 665. 
Budavari, T, et al., 2003, ApJ, 595, 59. 

Caldwell, R. R., Dave, R., Steinhardt, P. J., 1999, Phys. Rev. Lett., 80, 1582. 
Cooray, A., Hu, W., Huterer, D., Joffre, M., 2001, ApJ, 557, L7. 
Connolly, A., et al., 2002, ApJ, 579, 48. 

Eisenstein, D. J., in "Wide Field Multi-Object Spectroscopy", 2003, ASP Conference Proceedings, ed. A. Dey, 35. 
Eisenstein, D. J., Hu, W., 1998, ApJ, 496, 605. 

Feldman, H. A., Frieman, J. A., Fry, J. N., Scoccimarro, 2001, Phys. Rev. Lett., 86, 1434. 

Frieman, J. A., Gaztanaga, E., 1999, ApJ, 521, L83. 

Fry, J. N., 1994, Phys. Rev. Lett., 73, 215. 

Hu, W„ Phys. Rev. D., 2002, 65, 023003. 

Hu, W., Haiman, Z., Phys. Rev. D., 2003, 68, 063004. 

Hu, W., Jain, B., 2003, |astro-ph/0312395| 

Huterer, D., Phys. Rev. D., 2002, 65, 3001. 

Kosowsky, A., Turner, M. S., 1995, Phys. Rev. D, 52, 1739. 

Linder, E., 2003, Phys. Rev. D, 68, 083504. 

Linder, E., Jenkins, A., 2003, MNRAS, 346, 573. 

Ma, C. P., Caldwell, R. R., Bode, P., Wang, L., 1999, ApJ, 521, LI. 

Matsubara, T, Szalay, A., 2003, Phys. Rev. Lett., 90, 021302. 

Meiksin, A., White, M., Peacock, J. A., 1999, MNRAS, 304, 851. 

Mobasher, B., et al., 2004, ApJ, 600, L167. 

Perlmutter, S., et. al., 1999, ApJ, 517, 565. 

Pope, A., et. al., 2004, ApJ, 607, 655. 

Riess, A. G, et. al., 1998, AnJ, 116, 1009. 

Riess, A. G, et. al., 2001, ApJ, 560, 49. 

Seo, H. J., Eisenstein, D., 2003, ApJ, 598, 720. 

Seljak, U., Zaldarriaga, M., 1996, ApJ, 469, 437. 

Spergel, D. N., et al., 2003, ApJS, 148, 175. 

Takada, M„ Jain, B., 2004, MNRAS, 348, 897. 

Tegmark, M., et. al, 2004, Phys. Rev. D., 69, 103501. 

Tonry, J. L., et. al., 2003, ApJ, 594, 1. 

Verde, L., et al., 2002, MNRAS, 335, 432. 

Wang, L., Steinhardt, P. J., 1998, ApJ, 508, 483. 



© 2004 RAS, MNRAS OOO fnUTl 



